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ABSTRACT 

An electric field in the range of [o.3, 3.3] kV/mm 
is created normal to a thin-film FEP Teflon sample which 
accumulates potential of up to 8.8, 13.7 or 18.3 kV when 
exposed to an electron beam having energy of 10.0, 15.0 or 

20.0 kV, respectively. It is found that the secondary 
electron emission from the charged sample varies with field. 
The threshold voltage, at which the secondary electron 
emission coefficient o is unity, drops down from a low field 
value of 13.73 kV to a high field value of 13.11 kV for a 

15.0 kV beam- The maximum value of a measured at some 
primary beam voltage larger than the threshold falls off 
from 2.8 to 1.9. The beams are perpendicular to the 
sample's surface, i.e., parallel to the field. 

The established charge on the sample decays steeply 
near the edge of a metal-dielectric interface. A calcula- 
tion based on the potential shows a precipitous gradient 
of about 10 kV/mm being at the border. This gradient 
encourages the secondary electron emission as a charge 
transport mechanism by which stable charge distribution may 
be established near the interface. 

A computational technique has been developed that 
generates equipotential lines or contours and field vectors 
above a plane where potential is known. Furthermore, the 
utilization of conformal transformations allows the 
extension of the technique to configurations which map into 
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a plane. Nevertheless, it is assumed that the potential on 
the surface varies as a function of only one variable and 
that the equipotential lines or contours must start from the 
surface. Resulting lines or contours and vectors are two 
dimensional. 
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I. INTRODUCTION AND SURVEYS 


The foci of these experiments are charge accumulation 
and secondary-electron emission from a dielectric film in 
the region near a metal-dielectric interface which is 
irradiated by an electron beam. First of all, the charge 
(potential) distribution of that region was measured. 
Accordingly, an electric field pattern could be plotted from 
this information via electrostatic theory and computer pro- 
gramming. Then, the experiments examined the influence of 
field, both the computed field and an externally applied 
field, upon secondary-electron emission. 

Although we could not think of an effective way of 
measuring the secondary-electron emission near the inter- 
face, a qualitative attempt to deduce its behavior is 
made at the end of this report. .The deduction character- 
izes it in terms of the accumulated charge (potential) 
the calculated, immoderate field in that border, and in 
terms of the measured emissions at the central region where 
the undeflected beam is normal to the sample’s surface. 

Most of the experiments were with a FEP-Teflon sample, 
having thickness of 0.127 mm and one-side coated with silver 
inconel. This is one of the materials used to construct 
spacecraft. We also did some tests with 0.051-thick-Tef Ion 
sheet as well with aluminized Kapton which has properties 
similar to the Teflon. However, the thinner sheet of Teflon 
and the Kapton were found to be less suitable because of 
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leakages of electrons for high energy beams. In a vacuum 
chamber, having pressure of about lO"^ Torr, the sample is 
mounted under an aluminum aperture and exposed to an elec- 
tron beam with a kinetic energy of 10, 15, or 20 keV. The 
aperture not only defines the exposed area of the sample but 
also creates the metal-dielectric interface. A literature 
search of work done previously on the secondary— electron 
emission is briefly summarized in the succeeding paragraphs. 

The phenomenon of secondary-electron emission was 
discovered by Austin and Starke in 1902 [l], Austin and 
Starke were studying the reflection of cathode rays from 
metal surfaces and found that the metal target was able to 
emit a larger number of electrons than it was receiving. The 
striking electrons are called primary electrons and the 
liberated electrons are called secondary electrons. it is a 
complicated pheonomenon involving many different and intere- 
lated processes taking place at the target when its surface 
is irradiated by primary-beam kinetic electrons. There 
are three groups of secondary electrons emitted from the 
surface of a target: primary electrons reflected elastically, 

primary electrons reflected inelastically , and true secondary 
electrons. The energies of the secondary electrons released 
by a material vary from zero to the energy of the primary 
beam, and a typical energy distribution is shown in Fig. 

(I“l) foi* silver (solid) bombarded with electrons having 
energy slightly above 150 eV [2 J . N(E) is the relative 
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Figure (I-l) Energy distribution of secondary 
electrons. 

number of secondary electrons with energy E. There are 
three peaks in this distribution, denoted by a, b, and c, 
corresponding to the three types of emission: elastic, 

inelastic, and true secondary electrons, respectively. 

In application, people are often interested in the 
ratio of a number of primary electron over a number of 
secondary electron no matter- to which group the secondary 
electrons belong. This ratio, namely secondary electron 
emission coefficient o, may take any value from 0.5 for some 
metals to nearly 4.0 for some dielectrics. It -depends on the 
cleanness and temperature of the surface, on the angle of in- 
cident electrons, on. the work function of material, and of 
course, on the primary beam energy E . Cleanness, however, 

Jr 

does not cause any serious problem if the medium is kept at a 
sufficiently low pressure. It will be ignored here. 



4 


The fact that the secondary emission from metals is 
practically independent of temperature indicates that the 

r 

main process of energy loss is by interaction with the free 
electrons and not with lattice vibrations, which would 
certainly be temperature sensitive. In insulators inter- 
action of the secondaries with the lattice vibrations 
becomes significant, and this accounts for the reduction of 
secondary emission with increasing temperature observed in 
these materials. The coefficient a is high in insulators. If 
the primary beam is incident at an angle 9 to the normal to 
the surface, the secondary emission at low primary energy is 
very similar to the value for normal incidence; but at higher 
energy a increases as 9 is increased. The primary energy 
^prnax which a has its maximum value also increases as 9 
increases [2] . 

One might expect to find a high yield of secondary 

electrons from metals with a low work function, but this is 

not in accordance with the experimental results. McKay 

has plotted o as a function of the work function of 
max 

different metals and the result is just the contrary, 
namely a high work function corresponds to a high 
McKay observes rightly that the work function itself plays a 
relatively minor role in determining the secondary emission 
yield, but other physical properties such as the density 
are probably of more importance. When primary beam energy 
E is small the work function is apparently a governing 

hr 

factor [ 4 ]. 
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It is unnecessary to emphasize that the development 
of a theory of secondary-electron emission is no simple task. 
Various authors have made their contributions. Some theories 
are more or less phenomenological; others are based on a model, 
either classical, or the wave mechanical model of Bloch. The 
following paragraphs summarize. some of their investigations. 

According to Bruining [5], the liberation of secondary 
electrons occurs by the transfer of energy from the primary 
electrons to the electrons of the lattice. The behavior of 
the primary electron is determined by its energy loss as a 
function of penetration depth, its absorption, and scattering. 
The secondary electrons are scattered, and before reaching 
the surface a fraction will' be lost by absorption. Moreover, 
the work function has to be overcome before an emission is 
possible. He displayed an equation representing a universal 
curved based on theoretical considerations only, 


g 

g 

max 


1.85 F 


0. 92E 

— _ — E 

E 

pmax _ 


(i-i) 


where 


r 

_ 2 ^ 2 I 

F(r) == ® ^ dy , and r = i ~ 

o 


with, a and a depending on material. 

, It is striking that a universal curve, applicable to 
metals, can be found in actual fact as shown in Fig. (I~2) . 


ORIGINAL PAGE IS 
OF POOR QUALTIY 
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E /E 
p' pmax 

Figure (1-2) Universal o versus E_ curve: full 

line according to Baroody L 7 J, Eqn. 
(I-l) ; dotted line according to 
Jonker [ 6 1, Eqn. (1-2) ; the points 
represent measured data. 

We assume with Jonker [6] that the secondary . electrons are 
equally distributed in space at the point of origination. 
Then Eq. (I-l) must be modified to read: 


— ■ 


max 


pmax 


( 1 - 2 ) 


where 


1 V^/z 


f(r) = exp(-r^/z)z^ dy exp(y^), and r=Ey^ 


The smooth curve in Fig. (1-2) must be replaced by the 
broken curve, and thus is somewhat better accord with the 


observed data. 


til a 
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Second mention will be given to the theories of 
Kadyschevitz [ s] and later Baroody [ 9 ] , who used the 
Sommerfeld model and calculated the energy transfer to the 
conduction electrons along classical lines. Baroody did 
not follow his model consistently; he could have derived 
the relation between penetration depth and energy of the 
primary electron, but instead assumed Whiddington’ s law. 
Nevertheless Baroody made an important contribution to 
existing theories which has been discussed. 

Quantum— mechanical theories have been developed by 
Frohlich t lO] , Woldridge t 111 and Dekker and van Der Ziel 
[ 12 ] . These authors used Bloch model for metals and 
considered the collision of a single primary electron with 
a single metal electron. They paid special attention to 
the direction of the momentum transferred; thus in their 
model the emission of entirely free electron is impossible, 
since this would be contrary to the law of conservation of 
momentum. 

The most recent search on dielectric properties by 
Wall, Burke and Frederickson [ 13] states that although not 
all incident electron energies of interest have been covered 
for all materials, sound theoretical and semiempirical 
relationships have been developed that can be used to extend 
the available data. An example of such a relationship is 
the "universal secondary emission curve.” It is given by: 


max 


1-n 

X 

max 

g (X y 

'=’n' max 


E 


1-n 




-) 


(1-3) 


pmax 


Of POOH 
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where g (X) = t l-exp(-X)"J/x”“^ with X is the value of X 

max 

for g to be maximal. For a given material, X and n 
must be determined numerically to fit the available data. 
Most measured values of the secondary emission coefficient 
can be fit to the universal curve. In fact, if data is 
found that cannot be fit to the curve, there were probably 
errors made during the measurements. Fig. (1-3) shows 
secondary emission data for Teflon taken from Matskevich 
[ 14] fitted to the universal curve. The data was taken from 
a plot in the paper and deviations of some points from the 
curve are probably due as much to reading the plot as to 
experimental error. 


a 



Figure (1-3) o points are measured data for Teflon 
with . 21 , ®pmax~^ " ^ keV . Solid 

line is the curve yielded from Eqn,(I-3) 
with n=1.725. 



9 


For electron energies above about 0.5 keV, the 
following empirical relationship holds well: 

o = KE“® , (1-4) 

where K and m are constants. For most organics m is found 
to be about 0.725 and K depends on the specific material. 

The angular dependence of secondary emission follows the 
semiempirical relation: 

o„ = exp[ C(1 “COs6)3 , (1-5) 

“ o 

where 

0 = angle of incidence of electrons with respect 
to the surface normal 

= secondary emission coefficient at normal 
incidence 

Oq == secondary emission coefficient for electrons 
incident at angle 0. 

The constant C is determined empirically. For most polymers, 
C is approximately two. 

Because of its practical applications, secondary 
emission has long been a subject of exploration. As a 
result, a considerable volume of data exists covering many 
materials including metals and insulators. However, no 
theory has yet been found which satisfactorily covers all 
the observed phenomena. In as many books and publications 
as we could find, most of the authors deal with an 
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anode-cathode construction where the anode is a target, or 
with an arrangement where secondary electrons are liberated 
from the material itself by bombardment. The thing which 
we would like to point out is that we are studying the sub- 
ject of secondary-electron emission from a different view. 
The measurements are accomplished with a sample having an 
accumulated charge of electrons. They, therefore, include 
the effects of the accumulated electrons on the surface and 
the intrinsic electrons of the material. The second 
difference is the existence of an excessive field in the 
medium. Moreover, due to the unequally-distributed charge 
from the edge of the interface to the center of the charged 
sample the space above the sample possesses an inherently 
nonuniform field. A closely coordinated report [ by 
Robinson concentrates on the measurements of specimen 
charging and flashovers. 

Since apparatus and methods of measurement play a 
key role in attacking the goal, the descriptive Chapter II 
is for them. We then pay attention to the principal 
interests of this work, Chapter III and IV, field calcula- 
tions and secondary-electron emission in the presence of 
field. In Chapter III, we introduce a generalized-analytic- 
computational technique of equipotential and field mapping 
which are applicable to numerous electrostatic problems. 
Chapter IV is a mathematical formalism associated with the 
method of converting the outputs of the measurements, graphs 
on a recording chart, into a meaningful term, the 
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secondary-electron emission coefficient. It also displays 
plots of the coefficient versus beam energy and the 
coefficient versus estimated field. Chapter V is merely a 
review, and comments on important points in the preceding 
chapters. For clarifications, the appendix contains 
materials such as proofs of some expressions, tabulated 
formulae, and the computer program which was written for 
calculating equipotential contours. 



II. DESCRIPTION OF APPARATUS AND METHODS OF MEASUREMENT 


Apparatus should be appropriate for methods and vice 
versa to accomplish the objectives. Two major methods shall 
be applied, one being duplicated from Robinson's work [,16] , 
the other being a qualified modification of the first. 

Some pieces of apparatus were available on the market while 
some others were assembled in our laboratory to meet the 
methods' requirements. We tend to describe the performance 
of the apparatus rather than their dimensions which, not 
being crucial in the measurements, will be omitted for 
simplicity . 

A. Apparatus 

Apparatus is divided into three groups of equipment: 
the vacuum system itself, the apparatus inside, and the 
instruments outside of the vacuum. 

1, The vacuum system 

A chamber covered with a stainless steel bell jar 

having a diameter of 44.5 cm, and a height of 73.60 cm, is 

evacuated by a turbomolecular pump to a base pressure of 
0 

10 Torr. The. conventional vacuum system is equipped with 
an ionization gauge, a controllable leak, and assorted 
electrical and mechanical feedthroughs Fel, Fml, and Fm2 in 
the base plate as shown in Fig. (II-l) . Voltage feedthroughs 
Fe2 for the electron source are on the side of the bell jar. 

A glass window at the middle of the bell jar offers a 
view of the interior. 



tungste,n 
circle ^ 


aluminum 
^/Snap box 



Dual high 
(Voltage 
feed- 
f1 through 
Fe2 


Ti run. I 

nine pin feedthrough 
Fel 


rotary 

mechanical 

feedthrough 

Fm2 


Figure (II-l) Internal arrangement of the chamber 
has three constituents: an emitter, 

a regulator, and a target. 


ORIGINAL PAGE IS 
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2. Inside the bell jar 

There are three main components in the chamber : an 

emitter, a target and a regulator as shown in Fig. (Il-i) . 
The emitter is an electron source made of a 3,80 cm diameter 
circle of 10 mil tungsten in an aluminum snap box, as shown 
in Fig. (II-2) . Through the tungsten circle runs a 4-6 A 
current provided by the secondary coil of an isolation 
transformer where the primary voltage is adjustable up to 


120 volts. 



Figure (II-2) Structure of the emitter. 


The target essentially includes a specimen and three 
baffles. The specimen is mounted as shown in Fig. (II-l) 
on an eight-inch diameter aluminum disk which has a slot. 
Electrical contacts from the dielectric sample to the feed- 
throughs Fel underneath are established through this slot. 
The specimen defines the exposed area of the sample via an 
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aluminum aperture plate placed on top and creates a metal- 
dielectric interface. (Sample preparations shall be 
described in Section B.) Below the specimen, the baffles, 
surrounding the feedthroughs like curtains, prevent free 
electrons from striking the electrical terminals under the 
sample. 

The regulator consists of five elements: a collector, 

a screen, a shield, a baffle, and a shutter. 

a) The collector is an aluminum cylinder, having a 
radius of 51.5 mm and a height of 51.5 mm, being mounted on 
the aperture plate. This assembly is isolated from the 
grounding disk by the sample. Its function is to collect 
secondary electrons emitted from the sample and it is 
connected to an electrometer via the feedthrough Fel. The 
dimensions of the collector were chosen so that it would 
not disturb the charge accumulation on the sample. 

b) The screen with transparency of 85%, a 304-stain- 
less steel meshes made of 25.4 pm-diameter wire having 
0.2286 ram gap between each two of them, is hung above the 
sample and connected to the collector which is grounded 
through the electrometer to create a field between it and 
the charged sample. 

c) The shield is used with the collector. It is a 
three-layer triangle, shown in Fig. (II-3), attached to the 
rotary mechanical feedthrough Fml . 



attached to layer 1 

rotary feed- r> 

^through Fml / 



Figure (II-3) a. Top view of the shield 

b. Cross-section of the shield, the 
collector and the target (thick- 
ness of the sample is exaggerated) . 


Layer 1 is a grounded triangular sheet of aluminum 
on which a 2.6 ram inner diameter tube is positioned in such 
a way that only electrons whose trajectories are straight 
through the tube can reach the sample. Its sides are bent 
up to absorb the reflected electrons. Layer 2 is a 
dielectric sheet to separate layer 3 from the grounding 
layer 1. Layer 3 is a cap for the collector. It captures 
secondary electrons released from the sample which do not 
strike the collector. It is connected to the collector 
with a flexible wire. 

None of the dielectric surfaces of the shield are 
exposed to the beam or to secondary electrons. Otherwise, 
they would accumulate charge. The shield slides on the 
collector so that the tube projects a locus across the 



17 


sample. It can be slidden sufficiently far to completely 
expose the sample under the beam through the large hole on 
the baffle. 

d) The baffle is grounded and it prevents free 
electrons from striking the collector. 

e) The shutter is attached to the mechanical rotary 
feedthrough Fm2. it covers the tube, and the large hole on 
the baffle as well, therefore, eclipsing the sample when 
necessary. 

All the items are cleaned chemically with trichloro- 
lene and methanol before being put into the vacuum so that 
low pressure could be attained with a short pumping time. 
The sample, however, was not cleaned chemically, because 
chemicals might change its properties. 

3. Outside the bell jar 

Equipment outside the chamber consists of power 
supplies to the inputs, Ee2 , and meters to measure the out- 
puts, Fel. 

Two power supplies are needed. One provides an 
acceleration voltage to the snap box shown in Fig. (II-2) 
in the range of Eo.O, 22. o] Kilovolt d.c. and it biases the 
filament by the same amount. It is a partial supply for 
the emitter through feedthrough Fe2 . The other is a 
reference voltage in the range of the acceleration voltage. 
Reference and acceleration voltages are two terminals of a 
voltmeter where the difference is readable as low as 20.0 
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volt. The reference voltage is useful when we want to alter 
accelerating voltage by a small amount . 

Two electrometers designated Ml and M2 are used to 
register the responses of the dielectric sample and to 
amplify them in terms of coulombs, volts, or amperes. They 
are connected to the sample and the collector via Fel. A 
dual-channel recorder plots the amplified signals from Ml 
and M2 and its outputs are graphs on a recording chart. The 
electrometers and the recorder are calibrated adequately so 
that quantitative data can bo taken from the -graphs. 

B. Methods of measurement 

Two methods were employed. They are different in 
sample preparations, in equipment and in procedure. 

1. Preparation of samples 

Most of the work is with 5 mil FEP Teflon sheet which 
had been coated on one side with silver and then inconel. 

The material is mounted with the uncoated side exposed to an 
electron beam and the coating itself maintained at a virtual 
ground potential. The goal is to isolate an exposed region 
where we can do the measurements. The metal-film ground 
plane is cut into segments by an electrical discharge as 
illustrated in Fig. (II-4). Ref. [ isj is the origin of this 
fabrication. Each segment of the underlying ground plane 
was grounded by its own separate connecting wire. Measure- 
ments of charge (or current) were made by grounding a 
particular segment through an electrometer which measures a 
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Figure (II-4) a. A scheme to segment a metal coating 

of a dielectric sheet 
b. Top view of the specimen and the 
disk. 


charge (or induced current) while maintaining a virtual 
ground. Electrical connections to segments were constructed 
by bonding #32 copper wires to the segments with silver- 
filled epoxy. These wires lead to nearby ceramic insulated 
binding posts from which connections can be made easily. 

Two kinds of samples are fabricated to go with the two 
specific methods. Sample Sal has many small-rectangular 
segments whose top view was shown in Fig. (I 1-4) and whose 
cross-section is shown in Fig. (II-5) . Sample Sa2 has one 
circular segment for which the isolated ground region 
through the electrometer equals the area exposed under the 
beam (when the shield and the shutter are turned aside) . 
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Figure (I 1-5) Two samples shown with exaggerated 

thickness . 

a. cross-section of Sal. Charges on two 
segments are measured while other 
segments are grounded 

b. cross-section of the shield, the 
collector and the specimen with Sa2 
in action. 

2. Method Number 1 

The objective of this method [ig] is to gain a charge 
(voltage) distribution across a dielectric sample when it is 
being irradiated by a certain primary beam having a voltage 
of 6, 10, or 20 kV. We measure the charge on each 
isolated segment from the edge of the middle to the sample 
as well capacitance c^ of each. Then, the voltage is given 
by v^ = qjj^/CjL "til® center point of the ith segment. 

The equipment needed includes the sample Sal, the 
emitter, the target, two electrometers for measuring charges, 
the filament supply, and the acceleration supply. 
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Firstly, we check all the electrical connections, 
calibrate the strip-chart recorder, and connect two segments 
to the electrometers while other segments are grounded. Then 
we raise the filament voltage to produce the desired beam 
current. The charge measurements proceed as follows: 

a) Raise acceleration voltage to the desired value. 
Wait for a few seconds to charge the sample, but not so long 
that the accumulating electrons on the surface drift through 
the thin film. 

b) Register the maximal charge readings on both Ml 

and M2. 

c) Discharge the sample by decreasing the accelera- 
tion voltage to zero gradually. This step should be done 
slowly. Otherwise, there would be residual charges on the 
sample. 

d) Register the readings on both electrometers when 
acceleration voltage is zero. 

e) The accumulated charge for the segment connected 
to Ml under that beam voltage is the algebraic subtraction 
of the second reading (in d)) from the first reading (in b) ) 
on Ml. The charge on the segment connected to M2 is a 
similar subtraction of readings on M2. We may repeat this 
procedure a) - e) two or three times to be sure of the 
charge value obtained. Generally, it is nearly constant in 
the measurements for each segment. If not, we take the 
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average value. Charges for other segments are found by 
connecting the electrometers to corresponding positions in 
Pel and by repeating the procedure a) - e) . 

Capacitances are measured directly for each segment 
by placing a drop of aqueous salt solution, which acts as 

I 

an electrode, on the dielectric surface and applying 
potentials of 500 V and 1,000 V across the film. Corres- 
ponding charges q^ and are read on one electrometer. 
Usually, the capacitance for a particular segment is taken 
to be qj^/500 which is the same as q2/l,000. Capacitance 
measurements are done outside the vacuum. 

3. Method Number 2 

The objective of this method is to measure secondary- 
electron emissions from a sample, which has been already 
charged up with a charging beam, when it is irradiated by 
a probing beam acting as the primary beam. The approach is 
an effort to distinguish the two types of electrons, primary 
and secondary, via readings of currents which are induced 
on the sample and the collector. 

Sample- Sa2 replaces Sal. In addition to those things 
used in method Ir-the emitter, the target, the electrometers, 
the filament supply, and the acceleration supply': — w® need 
the regulator, the reference supply, the recorder, and for 
some specific measurements the fine screen which hung above 
the sample within certain distances from the sample. 
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As a preliminary, we check all the electrical 
connections calibrate the recorder and electrometers so that 
graphs on the recording charge are time varying plots of the 
electrometer readings. We then connect M2 to the sample Sa2 
(one segment only) , Ml to the collector and raise the fila- 
ment voltage to have the wanted beam current density. The 
secondary electron emission measurements are then made as 
follows: 

a) Open the shutter and the shield so that the 
sample is entirely exposed through the large hole on the 
baffle of the regulator. 

b) Increase acceleration voltage to the desired 
value of charging beam voltage. The sample is charged up 
freely as if the regulator were not there. 

c) Wait for a few seconds and then close the shutter 
and the shield in that order. The charged sample is no 
longer under the beam's irradiation. 

d) Decrease acceleration voltage to probing-beam 

voltage. This should be done quickly to avoid drifting of 

the accumulated electrons on the surface through the sample, 

-9 

e) Set Ml and M2 suitable scales (10 ampere full 
scale), and start the recorder. Steps d) and e) should be 
completed in less than 30. seconds, 

f) Open the shutter only, while the shield remains 
on top of the collector. Hence the unique path for the 
electrons that strike the sample is a straight line through 
the tube onto the shield. 
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g) Simultaneously with f ) , let the recorder draw 
instantaneous responses of the meters due to the induced 
currents on the sample (M2) and on the collector (Ml) . We 
repeat the procedure a) - g) with the probing beam voltage 
varying from a small level, say the primary voltage less 

2.000 volt, to the charging beam voltage in steps of 50.0 
volt measured against the reference voltage. 

Results from method 1 show us that the sample's 
surface voltage, namely the threshold voltage, is within 

3.000 volt of the primary beam voltage. Method 2 is based 
on the assumption that we have a reading on M2 (sample) 
whenever the secondary beam voltage is greater than the 
threshold voltage. Otherwise, the beam will be reflected 
by the negative surface potential and collected by the 
collector . 

We must accelerate the beam with two different levels 
in two unlike environments. The higher level beam, the 
charging beam, charges the sample up openly. The lower level 
beam, the probing (primary) beam, irradiates the charged 
sample with the narrowed beam passing through the tube. The 
tube is oriented to the center of the sample and it is small 
enough so that the spreading electrons do not strike the 
aperture or the collector. Therefore, as soon as they enter 
the tube, the primary electrons (probing beam electrons) 
must be in one of the two following instances: 1) They are 

reflected and collected if the probing beam voltage is 
smaller than the threshold voltage. There is no response on 
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electrometer M2 and a constant induced current on electro- 
meter Ml. 2) They strike the sample if the probing beam 
voltage is greater than or equal to the threshold voltage. 
Then secondary electrons are given off from the sample. 
Consequently, the responses occur on both electrometers, and 
therefore, two graphs versus time appear on the recording 
chart. The graph corresponding to the collector is the 
voltage induced by the emitted-secondary electrons. The 
graph corresponding to the sample is the induced voltage 
due to the secondary electrons minus primary electrons 
entering through the tube. Information contained in the 
graphs shall be further analyzed in Chapter IV, Section A. 

4'. Time constant measurement 

For calculation in Chapter IV, one needs to 

know the time constant of the circuit from the sample 

through the electrometer M2 to the graph on the recording 

chart. It is measurable by method 2 with the primary 

(probing) beam irradiating the sample for a short duration 

■t's (fraction of a second). The induced current ^ 2 ^^^ 

sample becomes a pulse. The graph becomes a pulse response 

of the low-pass filter as shown in Fig. (II-6) . On the 

graph we select one point (t^', v(t^)) to be the reference. 

The time constant is the difference t minus t^, with t 

satisfying v(t) = v(t )/e The time constant is measured 

r 

in the same environment where the secondary electron 
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Figure (II-6) a. A low pass filter with a pulse input 
b. One way to measure the time constant 
from the response: RC=t-t for t >t . 

emission happens. It varies from case to case, with the 
screen or without the screen and somewhat with the amount 
of charge accumulated on the sample's surface. 


OBIGINAL 
OF POOH QUAJJTYi 



III. FIELD MAPPING ABOVE A DIELECTRIC SURFACE 


After a few measurements with sample under circular 
apertures having radii of 2.5, 3.8 and 5.1 cm, the fact 
became apparent that the voltage gradient is unchanged near 
the edge as shown in Fig. (Ill-lb) . This leads us to the 
proposition that effects of walls other than those at A and C 
do not appear in the variation of potential along x-direc- 
tion, and that along z-direction illustrated in Fig. (Ill-la) 
the potential can be regarded as a constant. In other words, 
the surface’s potential changes with x only. 



(a) 


sample 


(b) 


Figure (III-l) a. The three-dimensional 

specimen 

b. An ordinary shape of voltage 
rapidly decays near the edges. 
For three different apertures, 
flat interval is longer for a 
larger radius sample. 
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A. Analytic Considerations 

In many practical situations, the boundary, either 
itself or after some complex transformations, is a plane 
with a potential specified as a function of- only one 
variable. Therefore, the region of interest which is the 
upper half space limited by the plane reduces a two- 
dimensional boundary value problem. The solution with- 
Dirichlet boundary conditions [it], which is proved in 
Appendix Al.l, .for an arbitrary position vector ^ becomes 

Ifr (v’;r-)ds’ 

s' 

where G(v;^') is a Green's function and the prime refers to 
the boundary. 

Suppose the boundary potential is given by 

4>(r') = V(x') , y' = 0 . 

Accordingly, the potential for y>0 is z-independent as shown 
in Figure (III-2a) . The Green's function is the potential 
at V = (x,y) due to a unit line charge at v ' = (x',,y'), 
while the x-z , plane is grounded as illustrated in Figure 
(III-2b) . Via the method of images, the total potential for 
y>0 is found to be 


G(v';T' ) = - Jin 


(x-x')^+(y-y')^ 
(x-x' )^+(y+y')^ 
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Figure (III-2) a. Dirichlet boundary conditions 

b. Method of images 

so that on the plane W = ^* • = (x',0), G(T;W) = 0 

the boundary condition- is satisfied [ isj . Subsequently, the 

normal derivative on S * is 


8 

an 


t-G(v;v’) 


(S’) 



■G(x,y;x’ ,y’) 


y’=o 


“4y 

/ 2 • 
(x-x') +y 


The final formulation of a potential at a point (x,y), y>0, 
is achieved: 


P(x,y) 



—00 


V(x’)y 

(x-x')^+y^ 


dx' 


(III-l) 
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The electric field is then derived from the potential; 


E(x 


00 

- ?7 


-oo 


V(x>) dx> 

C (x-x')2+y2]2 


•r2(x-x')y^ +[y^-(x-x')^]'e' 1. (III-2) 

X ^ ) 


where and are unit vectors , 

Let the voltage distribution on the boundary be 
piecewise linear, as shown in Figure (III-3) . Then the 



Figure (III-3) 


Piecewise linear voltage distri- 
bution for n points with constant 
wings. Vj and may be negative, 
zero, or positive. 
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linear functions of V(x’) are substituted into the 
expressions for P(x,y) and E(x,y) to yield the following; 


P(x.y) - |(Vi+V„) + (VjTj-V^T^) 
n-1 

+ (BjX+c.i)(T^^j-T^)+B.y(L,^j-L,)] 

i=l 

(III-3) 

E^(x.y) - (ViG^-V„G„) + 
n-1 

i=l 

E/x,y) - (ViPi-V„E„) + 
n-1 

r t (BjK+aj) 
i=l 

(III-4) 

where 

^ l[ (x!-x)^+y^] 

- (xi-x)D^ 
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Gi = y D^. 

, x! -X 

T. = — aJ’ctanf— i — ) 

IV y 

B. Discrete Numerical Approach 

Having established compact forms for P(x,y), E (x,y), 
and Ey(x,y) enables us to compute coordinates of points on a 
characterized equipotential contour and electric field which 
is normal to the contour at each point. For a given 
P(x,y) = P in the interval Cv^, ^i+1^ ’ always 

identify a point (x^^O) on the x-axis where 

P-«i 

^o = • (III-5) 

r 

Since Equations (III-3) and (III-4) do not hold when y=0, 
the subsequent point is assigned a small y value and has 
coordinates as follows: 


X 


1 


X 


o 


f 


yi 



i+l 


—IT* 


200 


(III-6) 


One calculates the field at this point and then traces along 
it, as illustrated in Figure (III-4) , to obtain 
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(new) 


y^(new) = 


E 


+ D 


+ D 


X 


2 2 


E„ 


E^+E^ 
X y 


(III-7) 


where D = P(Xj^,y^)-p controls the size of the step- 



Figure (II 1-4) A combination of steps, perpen- 
dicular and parallel to the 
field, that converge to a new 
point on the P contour (shown 
dashed) . 


The new pair (x^,yj^) is substituted into Eqn. (III-2) 
to check whether the difference between P and is 

less than some criterion of convergence of not, i.e.. 



(III-8) 
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where PSTEP is some constant depending upon the range of P, 
Steps parallel to the field are repeated as necessary. 

The next point is found by making a step perpendicular 
to the field to obtain 


X, 


I ^ 

- A I PSTEP — ./ ■ 

^ E^+E^ 

X y 


yj - yj + A|PSTEP| 

V®y 


(III-9) 


where -1 < A < 1 is used to select the size of the step; its 
sign is opposite to that of Field ECxgjyg) is then 

computed and Eqn. (III-7) is applied iteratively until the 
criterion, Eqn. (III-8) , is satisfied. This process is 
repeated over and over until we gain the desired number of 
points or. until we enter a region where E is very small. 


C. Experimental Results 

1. A sample under a circular aperture 
For a 13 mm radius sample of FEP Teflon with a 20 kV 
electron beam, the voltage distribution measured on the 
sample’s surface is given by Table (III-l) where u is the 
distance from the left edge. The geometry shown in Figure 
(III-5) can be converted to a plane with only one conformal 
mapping if the right hand edge is assumed to be at infinity. 
Letting z = z+iy and w = u+jv, one finds that the appropriate 
transformation [ 19] is 
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W 


= -i [ s+in (s+z) ] , 


4 




(III-IO) 


Also, Eqn. (III-IO) is employed with real z and w to find 
the points on the z-plane (y=0) corresponding to the data 
shown in Table (III-l) . 

A brief review of Eqn. (III-3) and (III-4) suggests 
that we linearly interpolate the data, and that the set of 
coefficients is given by 


V -V 

r> i + 1 i n 

’ “i - 


i = 1, . . . ,n-l 


(III-ll) 


The two wings are given as 


<! 

II 

<J 

II 

o 

for 

X < X 

V = = -17.0 

for 

X > X. 


1.0 

27.43 . 


The method described in Section B was carried out to find 
the contours from -2 to -16 kV in steps of -2 kV- The data 
points describing the contours were converted back to the 
original w-plane via Eqn. (III-IO). 

For the field components, the formalism of conversion 
is somewhat different. Since they are analytic functions of 
(x,y) we invoke Cauchy-Riemann conditions to relate field 
components E^, on the original plane to the field 
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Figure (III-6) Equipotent ial contours near a charged FEP Teflon 

surface in a 20 kV beam. 


CO 

<1 



38 


components E^, Ey on the transformed plane. It is shown in 
Appendix A1.2 that 




where 

f » 


dz 

dw 



) 


(III-12) 


is obtained from Eqn . (III-IO) . This sequence of Eqns. 
(III-l) to (III-12) and Table (III-l) were implemented. with 
a computer as illustrated in Appendix A3. The outputs are 
plotted in Fig. (III-6) . 

2. Sample under two rectangular apertures 
separated by a hemicy Under 

Similarly to the sample under a circular aperture, we 
first extend the outer boundaries to infinity, as shown in 
Fig. (III-7) (providing that the hemicylinder ' s radius is 
small compared to the sample dimensions) and then straighten 
the hemicylinder out via the transformation [ 19 ] 

z = w + . (III-13) 


For the field components, we use 


f • ■= 



1 


(III-14) 
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V 



Figure (III-7) Transformation of the system 

with a hemicylindr ical ground 
strip. 

Again, the sequence of Eqns. (III-l) to (III-12) , where 
(III-13) and (III-14) take the place of (III-IO) and 
(III-12) , respectively, were carried out to find the equi- 
potentials and field. Table (III-2) displays the data used 
for this example and Fig, (III-8) shows a plot of the out- 
puts . 



y 



16 kV 


Figure (III-8) Equipotential contours, near an FEP Teflon surface in 

a 15 kV beam. The grounded hemicylinder divides the 
surface symmetrically. Other points of negative x 
are symmetric about y-axis with points of positive x. 
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Table OCII-2) Data points representing surface 

potential; The specimen is symmetric 
and so is the voltage on its surface. 


Potential (kV) 

Coordinate, 

Transformed 

Coordinate 

-17.4 

-7 . 144 

-7.284 

-17.2 

-5.661 

-5.837 

-16.1 

-3.754 

-4.021 

“14.8 

-2.907 

-3.2511 

-11.6 

-1.636 

-2 , 247 

0.0 

-1.000 

-2.000 

0.0 

1.000 

2.000 

-11.6 

1.636 

2.247 

-14,8 

2.907 

3.251 

-16.1 

3.754 

5.837 

-17.2 

5.661 

7.284 

-17.4 

7.144 




iv. INFLUENCE OP FIELD UPON SECONDARY ELECTRON EMISSION 


In the preceeding chapter, we portrayed equipotentials 
^jjd fields calculated from, the voltage distribution on the 
gg,niple's surface which had been measured via method 1. In 

chapter, we shall calculate secondary electron emission 
coefficients for which data were taken via method 2. The 
^jj-terpretation involves solving an inhomogeneous differential 
ea^ation where the right-hand side is the sample’s substrate 
current source. The topic of secondary electron emission 
field is examined in the last section. 

A . Analysis 

Based on the continuity relationship between current 
ajid charge, the following arguments deal with currents (or 
induced by charge (electron) fluctuations. What 
could see are outputs of two parallel RC filters on a 
j^ccording chart as illustrated in Fig. (IV-1) . There are 
always two curves Vj^(t) and v^(t) associated with each 
measurement. 

The response v^(t) of the circuit is produced 

the current ij^Ci) flowing from the collector being 
connected to the electrometer Ml. When secondary electron 
emission reaches its steady state, the number of primary 
gXectrons striking the sample equals the number of secondary 
electrons released from the sample. The current i^(t) 
liecomes a constant and then v^(t) = Ril^ is the voltage 

created by primary electrons. The response ^ 2 (t) of the 




steady state 


(b) 


a. Schematic of recording 
circuit. 

b. Typical outputs of two 
RC circuits on the 
chart . 
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^2^2 to the induced by the 

number of electrons from the sample surface, 

secondary electrons less ^ 3 ^tniary electrons. it is the 
solution of the first ord^^ differential equation 

l/Rg VgCt) + Cg dv2<^t;)/dt = i2<t) 


After several attempts, arrived at a physically rear 

sonable assumption for ig » that is 12 (f) = l 2 e~^'^^ , 

t>0. In a sense, the emi^^^*^“ from the surface is greatest 
at the instance t = 0 when first irradiate the surface 

which was .already charged- There are two unknowns, i and t 
with Ig representing the nia-ximum difference between numbers 
of electrons arriving and Xe^^ving at the sample's surface. 

We solve the differential equation with the condition that 
there is no stored energy ioiiislly in the capacitor, 

V 2 (o) = 0. The solution form 


V2(t) 


^ 2 ^ 2 ^ 

TE^cpFr 


-t/R2C2 

e 



(lV-1) 


where R 2 C 2 measurable a^d varies with cases. 

Even though the fofcius function is strongest at 


t = 0, the response is maJcinium at t^>0 due to a delay of 
the R2^2 i.e,, 
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d 

dt 


V2(t) 


'm 


I9R9T 


-e 

_ ^2^2 


+ 



= 0 , 


(lV-2) 


where and V 2 Cfcj^) are read from the chart. Substituting 
tm into Eq. (IV-2) , we gain the root t by using a bisection- 
iterative method.- Using t, t^^, and. V 2 (t^^) , we compute 
^2^2 (IV-1) . With R^I-j^ explained in the foregoing 

paragraph, the secondary electron— emission coefficient o is 
the ratio 

^ ^ number of secondary e lectrons _ ^l’’"^2 
number of primary electrons * 


We can choose the scales of both electrometers so that 
^2 ^ ^1‘ Consequently, a can be written in terms of chart 
voltages as follows: 


a 


^2 ^^l'*'^2^ 


1 + 


^2^2 

Rlll 


(IV-3) 


where ^2^2 computed from Eq. (IV-1) and Eq. (IV-2) , and 
Rflf is the steady state voltage due to primary electrons 
alone read on the plot of Vj^(t) , 


B. Analyzed Data 

On each two-graph output shown in Fig. (IV-1), we 

measure three factors t , v„{t ) and RtI, to be used in 

m’ 2' m 11 

OEIGINAI- 
OEPOOR 
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-three equations (IV-I) to (IV-3) . The coefficients versus 

primary beam voltage, Figures (IV-2) to (rv-6), were computedV^ 

for a 5 mil FEP Teflon sample under a circular aperture hav- ^ 

ing radius of 16 mm with charging beam voltages 10, 15 and 

20 kV. The difference between the two beam voltages is E » 

P 

Looking through Figures (IV~2) to (IV~6) we conceive 
that the curves are rapidly increasing with primary (probing) 
beam energies larger than the threshold voltage by a small 
amount. Within a 200—volt increment, they reach peak values. 
Then they smoothly decay as the primary beam gets larger. 

The threshold voltages are always at least 1.2 kV less than 
the charging beam voltage. If we let the charging beam 
voltage be the origin, then thresholds of beams of 20,0, 

15.0, and 10.0 kV appear in this order with respect to the 
probing beam. The peak values of the of the coefficient are 
about two when the screen is low and about 2.8 when the 
screen is 23 mm high above the sample surface. Note that 
in all cases, o is a very sensitive function of secondary 
beam voltage in the interval around c 

max 


Q, Accuracy 

Looking at Fig. (IV-1) closely, we realize v„(t ) 

2 ' m 

can be chosen more precisely than t^^. Since there is a 
range where we can select a value for t^, one question 
arises concerning the sensitivity of a with respect to the 
choices of and VgCt^^^), as well with respect to the 
measured ^ 2^2 inexactness of the assumed exponential 

form for igCl) la analysis section. 
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It was numerically found that 


A A 

^^0.8 ^ ^ 


AV2(t ) 
2 ' 


AO ^2^2 ■ 


where the normal values of t and R„C« are respectively 

in A ^ 


= 0.29i0.025 sec and = 0.155±0.01 sec 


Suppose the function i 2 (t) produces a 5% error on o, the 
largest possible deviation of a is then 


A 

— < 0.8 -:p- 

" - 'in 


m ^^2^2 

-+ 0.02 + 0.6 -^s=yA + 0.05 = 0.17768. 

*^ 2'"2 


We would claim that accuracy of this interpretation is 
within 18%. 


D. Field Influence on the Secondary Electron Emission 

As a matter of fact, the field had influenced the 
secondary-electron-emission coefficient computed in Section 
B. We ought to recognize two kinds of field being present 
inside the collector where the primary electrons and 
secondary electrons were captured: The accumulated field 

is the field due to charge accumulation on the dielectric 
surface. The enforced field is the field due to the 






screen's presence. 
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At the center of the sample toward where the beam was 
deliberately directed, the field possesses only a y-component 
because of symmetry of the boundaries. When the screen is 
near the sample, the enforced field predominates over the 
accumulated field and the total field is mostly uniform. 

When the screen is high above the sample, there is a fringing 
effect on the enforced field and the accumulated field was 
about equal to the total field. There also exists an elec- 
trostatic force which pulls the screen down, i.e., the 
distance between the screen and the sample is hard to con- 
trol. Various contributions make it difficult to obtain 
exact values of the field. Nevertheless, approximations of 
the field ranges with respect to the screen's position and 
surface (threshold) voltages tender us an idea of how the 
characteristics of secondary electron emission vary with 
field, as illustrated in Table (IV-1) and in Figures (IV-7) 
to (IV-9) . 

Fig. (IV-7) shows the differences of the a curves, 
with and without the field. As a matter of fact, there are 
scattered points due to errors made during the measurements 
and the computation. Scattering is greater when the field 
is stronger and also the secondary emission peaks are 
broader. Clearly from Fig. (IV-8) , the voltage interval 9 
between unity crossing points becomes larger with a higher 
field existing above the surface. When the screen touched 
the sample, we observed that surface potential could not reach 
its equilibrium as usual even for a 10 kV beam. The surface 



Table (IV-1) Estimated range of E„ = Eg+E^ yielded 
from the compound of^accumulated field 
Eg^ (in Chapter III) and the enforced 
field Ee (surface potential)/(screen ' s 
height) , 


O O 

o i 

s § 


Pr imary 

Screen's distance 

Threshold 

Approximated 


probing 

from the sample 

voltage 

field range 

Note 

beam (kV) 

surface (mm) 

(kV) 

Ey (kV/mm) 



4.0 ± 1.0 
10.0 ± 2.0 

18.15 

1.82 ± 0^20 


20 

15.0 ± 2.0 

18.25 

1.22 + 0.25 



23.0 ± 2.0 

62.0 ± 0.0 

18.27 

0.80 ± 0.30 




4.0 

+ 

1.0 

13.11 

3.28 

± 

0.15 

14.5 kV beam 


10.0 


2.0 

13.26 

1.33 

+ 

0.20 


15 

15.0 

± 

2.0 

13.40 

0.90 

± 

0.25 



23.0 

± 

2.0 

13.47 

0.59 

± 

0.30 



62.0 

± 

0.6 

13.73 

0.22 


0.45 

without the 









screen 


4.0 


1.0 

8.18 

2.05 


0.15 



10.0 

± 

2.0 






10 

15.0 

1 

2.0 

8.43 

0.56 

+ 

0.20 



23.0 

T 

i 

2.0 

8.50 

0.37 

i 

0.25 



62.0 

0.0 

8,75 

0.15 

+ 

0.45 

without the 









screen 


Oi 



Figure (IV-7) Data for 15 kV beam only with different positions 

of the screen: ><. without the screen, 23 ram, 

A 15 mm, oio mm, and o 4 ram from the sample’s 
surface . 






Figure (IV-8) Plots of voltage interval 6 versus magnitude 

of estimated field for xlO, o 15, and & 
20 kV beams . 
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accumulated charge to some level and then suddenly dis- 
charged. Is this due to a very strong field (zero distance 
of the grounding screen from the charge surface) or an 
extremely low threshold voltage? We doubt that discharging 
is caused by secondary electron emission. It is more likely 
another phenomena which is beyound our scope. For a 
discussion about discharging (flashover) refer to the report 
of Robinson [.15j , 

It is logical to think that under the influence of a 
stronger field the number of secondary electrons liberated 
from the charged- surface by bombardment would be larger. 

This is, however, not true. Rather, a is seemingly a 
monotonic-decreasing function of field and the sharpness 
around a on the a-curve is reduced by a high field as 
shown in Fig. (IV-7) and Fig.. (IV-9) . 



V. SUMMARY AND DISCUSSION 


As an analogy, the whole system can be depicted as a 
two-port network shown in Fig. (V-1) , where the "black box" 
is the vacuum containing a dielectric sample, an electron 
source, and some apparatus to support the measurements 
realized outside it. The inputs, filament voltage and 
acceleration voltage, supply power for the electron gun. 

The outputs, accumulated charges (or voltages) or induced 
currents (or voltages) , are read on the electrometers or on 
a recording chart. Inside the vacuum, when the sample is 
being irradiated by an electron beam, at least two processes 
take place: charge accumulation on the sample's surface and 

secondary emission from the sample's surface. They both were 
qualitatively and quantitatively determined. 
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Figure (V-1) A schematically equivalent portrayal 

for the experiment. The parenthesized 
numbers are ranges in which we 
operated. 
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Two objectives were achieved via two methods. Method 
1 [ 16] with many small-rectangular segments yields an 
accumulated charge (or voltage) distribution across the 
sample. From this distribution we solved a boundary-value 
problem. With the aid of a computer, numerical values of 
potential and field above the sample were found. As a 
continuation of the work, method 2 with a sample having a 
single-circular segment was used to measure the secondary 
electron emission from the sample's surface which was 
exposed to these fields and the fields enforced by the 
screen's presence. Information about the emission is 
conveyed from the collector and the underlying plane sample 
to the graphs on the recording chart by the induced currents 
on them. The interpretation requires mathematical manipula- 
tions and computer programming. 

The apparatus employed in method 2 is that employed 
in method 1 with some modifications. However, the rectang- 
ularly segmented sample which is hard to fabricate becomes 
unnecessary for the secondary-electron emission measurements. 
In contrast, method 2 is inapplicable for charge measurements 
near the edges of the sample, because the narrowed beam would 
be largely deflected when it is directed toward them. 

In order to see what the field patterns and equi- 
potential contours look like in the vicinity of dielectric 
samples where charge accumulation appears, two different 
geometries were investigated. One was a sample under two 
rectangular apertures separated by a hemicylinder . In 
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accordance with experimental results, and for mathematical 
convenience, voltage distribution on either sample's surface 
could be treated as a function of one variable. 

The potential everywhere in the upper-half space 
limited by the specimen satisfies Laplace's equation with 
Dirichlet boundary conditions. Applications of Green's 
theorem and Green's function simplify the solution 
significantly. The consequent potential and field are 
integral expressions of two variables. These analytic 
expressions combined with a so-called "trace the field" 
concept are powerful for numerical purposes. Fundamentally, 
tracing the field is performing a sequence of steps in the 
correct directions, parallel or perpendicular to the field 
at each point, under the control of certain convergence 
standards . 

A computer subroutine (in Appendix A3), which imple- 
ments the Green's solution for the Dirichlet boundary value 
problem and the concept of tracing the field, was written. 

It linearly interpolates a discrete set of boundary-potential 
points on the x-axis of a plane where potential is unchanging 
along the other axis. Then it produces equipotential 
contours and field patterns above that plane. The equi- 
potentials and fields are functions of x and of the vertical 
variable in the three-dimensional space. Conformal trans- 
formations enable this subroutine to spread its usefulness 
to any electrostatic problem whose boundaries can be mapped 
into a plane. The two restrictions are that the potential 
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varies with only one variable on the plane and the contours 
must start from the plane. 

In the measurement of secondary-electron emission 
from the .dielectric sample by a probing beam when the 
sample had been already charged up by a charging beam, two 
phases were involved. The first phase was a strategy of 
distinguishing between primary electrons striking the sample 
and secondary electrons reaching the collector. Then the 
electron induced currents, amplified with the electrometers, 
were registered simultaneously on a recording chart. The 
second phase was merely the solving of an inverse-filtering 
problem. The forcing function to be specified is the 
current signal induced on the collector. The circuit is a 
simple parallel RC filter. The secondary-electron emission 
coefficient is the ratio of the maximal-secondary-electron- 
induced current to the constant-primary-electron-induced 
currents. In brief, the general behavior of the coefficient 
with respect to beam voltage in this experiment is analogous 
to the result accomplished by many other authors as shown 
in Fig. (V-2) . More rigorously, as we should perceive, some 
factors associated with the secondary-electron emission do 
vary with an electric field. 



Figure (V-2) An ordinary curve for a with 
factors: Omax? 6) 
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With the screen's presence, it is obvious that the 
closer to the sample the screen is the stronger the field 
existing on the sample's surface becomes. Even though we 
were unable to ascertain precise values for the fields, 
comparisons shown in Fig. (IV-2)-Fig. (IV-6) among the 
threshold voltage intervals 9, among the peak values o , 
and among the deviations 6 do confirm the following state- 
ments ; 


a) The threshold voltage (surface potential) , at 
which the number of primary electrons equal the number of 
secondary electrons, is smaller for a larger field. There- 
fore, the threshold interval 9, the voltage distance 
between two unit-a points, is larger for a greater field . 

b) Secondary-electron emission coefficient a is a 
sharply increasing function with respect to primary voltage 
when it is slightly greater than the threshold voltage. We, 
therefore, may find more precise values of a with a finer 

1U2.X 

division in the primary beam. Anyway, the values of a 

niRX 

become low for a high field except at three data points where 
we presume that errors were made during the measurement and 
the calculation. 

c) Without the field, the coefficient a decreases 
70% with a 50-volt shift to the left and 40% with a 50-volt 
shift to the right of the primary voltage from the voltage 
corresponding to peak value a of the a curve. In the 

IUaX 

presence of the strong field these percentages approach 20% 
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and 10%, respectively. In other words, if we call deviation 
6 the voltage interval between two points on the o-curve 
at which a decreases a half with respect to unity from its 
maximal value, then the deviation 6 is larger with an 
influence of a stronger field. 

Results from. the computation in Chapter III reveal 
that the x-component of the field does not change much 
relatively compared to the y-component from one point to 
another across the sample, and that the x-component 
intensively dominates the y-component at field points near 
j the edge. This revelation is merged with conclusions in a), 
b) and c) to deduce an explanation why the surface (thres- 
hold) potential dies off abruptly near the metal-dielectric 
interface. If we assume that the tangential field has an 
influence on secondary electron emission similar to the 
normal field, then we may conclude that because the steep 
gradient exists in the border, the threshold voltage must 
be low. The accumulated electrons are always motivated by 
this field to move to the edge. Therefore, the secondary 
electron emission (obviously, with the coefficient being 
reduced but greater than unity) is readily excited even at 
low surface potentials by the charging beam. Consequently, 
the gradient is maintained by a balance of the charge trans- 
fer process of secondary-electron emission with the primary 
electron beam. This discussion is consistent with the 
pattern of charge accumulation on a dielectric sheet which 
was firstly seen in the experiment, but was not explained. 



APPENDIX A 


A1 Proofs of Expression 

1. Solution of Diriclilet problem with Green’s 
function 

Starting from the Green's theorem: 

=y[c}> |i _ ^ Jds , 

V s'" 

a person is able to convert the Poisson differential 
2 

equation, y 4“ = -4-Trp , into an integral equation if he 
chooses a particular i|/j namely Green’s function 
G(v,v’) = l/|v;v’| where v is the observation point and 
is the integration variable. Further, we put 4*=^, a scalar 
potential, and make use of 7 G(v;v’) = -4ir6(v-v’), so that 
the theorem becomes. 


-4tt^(V'’ ) 6 ) 

V 

+ 4ttp (W)G(^;^' ) ]'dv' = ^ )-g~ 7 -G(^;T' ) 

S 

- G(^;7^')-g~(7;7'')]ds' . 


If the point v lies outside the surface (S) , =0. If 

V lies within the volume (V), then the equation becomes: 
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p ) dv ' 

V 




- G(v; v’)^^(v;v’)l ds 


Two remarks are in order about the resulting expression 
. First, if the surface (S) is at infinity and the 
electric field on (S) falls off faster l/|^wj, then the 
surface integral vanishes. Second, for a charge free 
volume the potential anywhere inside the volume is 
expressed in terms of the surface potential and its normal 
derivative. 

The solution of the Poisson or Laplace equation in a 

volume (V) with either Dirichlet or Neumann boundary 

conditions on the boundary surface (S) can be simplified 

by choosing an appropriate = l/j( )j+F , 

2 _<■ 

where ^ F(v;v') = 0, to eliminate one or the other of the 
two surface integrals. For Dirichlet boundary conditions 
we demand that GC^;^) be zero when W on (S) . Then the 
solution is: 




V 


4ir 


/ 




We had a charge free volume such that the first term 
vanishes . 


obigwal pag^ 
OF POOR 



67 


2. Complex transformations for field components 
Let the variable z = x+jy relate to w = u+jv by 
z=f(w), where f is the transformation from the z-plane to 
the w-plane. Let the function y be given by y = 
where 4 is potential and i{/ is flux. Since tf’ and ^ are 
analytic function of (x,y) and therefore, of (u,v), the 
Cauchy-Riemann conditions permit the following manipula- 
tions : 


E. 


X 


and 


E 


u 


94> 

9\j/ 

; E . - 



9x 

y ^ 1 

X y 

X 

9x 


9'}’] 

9;}/ 

E == ~ I 

- ^ 

9u 

9v 

’ V 9v 

9u 


V 

u I 
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Let's consider 


9y 

Ju 


- ~ “E„+jE„ , 


U V 


and 


9y 

9x 


9<i> . 94 / ^ 


Then 9y/9u can be written in terms of E^, Ey and the 
derivative of the transformation as follows: 


_! r 9y ay aiyi 

2 .9x'9u 9y‘9u 




9u 
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Equating the real parts and imaginary parts of the two 
equations of dy/du, we construct the expression as used. 


A2 Tabulated Formulae 
1. Integrals t 20l 


/ 


xdx In. ,2. 
— T = 2b ^«(a+bx ) 
a+bx 


I 


xdx 


-1 


(a+bx^) 2bm(a+bx^)™ 


for m>l 


X dx 


-X 


4 * 


(a+bx^)®'^^ 2mb(a+bx^)"^ (a+bx^)”^ 


/■ 


dx 


dx 


2,22 

a+bx 


1 tan'l M 
ab a 


/■ 


xdx 


(a+bx) 


Un(a+bx) + 


These were used to evaluate the integrals of P(x,y) E (x,y) 

X 

and Ey(x,y) in Chapter III. 

2. Green’s function [4] 

For z-independent three-dimensional space, the Green's 
function has the form: 


G(^;W) - -2 J^n ( j j V*’ I ) , where 




V = (x,y) and v' = (x',y') . 
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3. Conformal mapping [l9J 

a) w = f(z) = i[s+ln(s+z)] , s = \jz^-l 


b) w = f(z) = z+l/z . 


A3 Computer Programs — Limitations of the 
Generalized Computer Program 


Generally, this program is able to find the equi- 
potential contours and field components of any configura- 
tion which can be converted into a two-dimensional space 
where the boundary voltage on the horizontal x-axis is 
well-specified in the form of Fig. (III-3) . Nevertheless, 
in order to use it efficiently one should learn several 
features of its structure. 


The two main programs were used to call subroutine 
Green (below) in Chapter III are; 


1. Main Program 1 


COMMON /ORI/ C(30), V(30) , BE(30) , AL(30), PY, NN 
COMMON /ARR/ XG(10 , 50) , YG (10 , 50) , FX (10 , 50) ,FY(10 , 
50) ,PP(10) ,KK(10) 

READ (5,1) N, NPOTEN, NPOINT, PMIN, PSTEP,ADJ 

1 FORMAT (3I5,4F10.0) 

PY = 3,1415926 
NN'= N - 1 

DO 10 I = 1, N 

READ (5,2) C00R,V(I) 

2 FORMAT (2F20.0) 

C(I) = TRNSFM(PY,C00R) 

WRITE(6,3) V(I) ,C00R,C(I) 

3 FORMAT (IH , 3 (20X, E13 . 5) ) 

10 CONTINUE 
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CALL GREEN (M, NPOTEN, NPOINT , PMIN, PSTEP, ADJ) 

DO 101 I = 1, M 

WRITE (6,4) PP(I) 

4 FORMAT (///,20X, ’P- ' , 1 (E13 . 5 , 20X) ///20X, 

’ABSS’ ,26X, ’ORDI' ,26X, ’EX' ,28X, 'EY'///) 

IF ( KK(I) .LT. 0 ) KK(I) = -KK(I) 

JJ = KK(I) 

DO 102 J = 1, JJ 

CALL INVRSl (PY,XG(I, J) ,YG(I, J) ,X,Y) 

IF ( J -EQ. 1 ) GO TO 103 

CALL INVRS2 (PY,X, Y, FX(I , J) , FY(I , J) , EX, EY) 
WRITE (6,6) X, Y, EX, EY 

6 FORMAT (IH , 4 (17X, E13 . 5) ) 

GO TO 102 

103 WRITE (6,6) X, Y 

102 CONTINUE 

IF ( KK(I) .LT. 0 ) PRINT 7 

7 FORMAT (/20X, 'SUBPROGRAM SCAN COULD NOT FIND THE 
1, SUCCESSIVE POINT’) 

101 CONTINUE 
STOP 
END 

FUNCTION TRNSFM(PY,W) 

Z = 1.2 

IF ( W .EQ. 0.0 ) Z = 1.0 
DO 30 I = 1, 30 

F = (1/PY) * (SQRT(Z**2-1)+ ALOG(Z + SQRT(Z**2 - 1))) 

1, - W 

IF ( ABS(F) .LE. l.OE-7 ) GO TO 101 
G = (1/PY) * ( (Z+1)/SQRT(Z**2-1) ) 

ZIT = Z - F / G 

IF ( ABS(ZIT-Z) .LE. l.OE-7 ) GO TO 100 
Z = ZIT 

IF ( ZIT .LE. 1.0 ) ZIT = ( (Z-l)/2 + Z ) 

30 CONTINUE 

100 TRNSFM = ZIT 
RETURN 

101 TRNSFM = Z 
RETURN 
END 

SUBROUTINE INVRSl (PY,ABSS, ORDI ,X,Y) 

COMPLEX U, Z, W 
Z = CMPLX( ABSS, ORDI) 

U = CSQRT( CMPLX((ABSS+*2-0RDI**2-l) , (2+ABSS*ORDl)) ) 

IF ( ABSS .LT. 0.0 ) U U 

W = (1/PY) * ( U + CLOG(Z+U) ) 

X = REAL(W) 

Y = AIMAG(W) 

RETURN 

END 


oeiginm-wg^ 
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SUBROUTINE INVRS2 (PY, X, Y,FX, FY, EX, EY) 

COMPLEX Z 

Z = CSQRT( CMPLX((X*+2+Y**2-l) ,-2*Y) ) / SQRT( (X-1) 
1 =k=i<2+y**2) 

G = REAL(Z) / PY 
H = AIMAG(Z) / PY 
EX = FX*G + FY*H 
EY = -FX*G + FY*G 
RETURN 
END 


2. Main Program 2 


COMMON /ORI/ C(30), V(30) , BE(30) , AL(30) , PY, NN 
COMMON /ARR/ XG(10 , 50) , YG(10 , 50) , FX(10 , 50)FY(10 , 50) 
1,PP(10) ,KK(10) 

W(Z) = Z + 1 / Z 

READ (5,1) N, NPOTEN, NPOINT, PMIN, PSTEP,ADJ 

1 FORMAT (3I5,4F10.0) 

WRITE (6, 8) 

‘8 FORMAT(///,21X, ’ABCISSA' ,15X, 'MODIFIED ABCISSA’,25X 
1, 'POTENTIAL' ,///) 

NN = N - 1 
DO 10 I = 1, N 

READ (5,2) COOR,V(I) 

2 FORMAT (2F20.0) 

C(I) = W(COOR) 

WRITE(6,3) COOR,C(I) ,V(I) 

3 FORMAT (IH , 3 (20X, E13 . 5) ) 

10 CONTINUE 

CALL GREEN ( M, NPOTEN, NPOINT, PMIN, PSTEP, AD J) 

IF ( M .EQ. 0 ) STOP 
DO 101 I = 1, M 

TOITE (6,4) PP(I) 

4 FORMAT (///,20X,'P= ' , 1(E13, 5,20X)///20X, 

1, 'ABSS' ,26X, 'ORDI ' ,26X, 'EX' ,28X, 'EY'///) 

JJ = KK(I) 

IF(KK(I) .LT. 0 ) JJ = -JJ 
DO 102 J = 1, JJ 
IF(J-EQ.l) GO TO 103 

CALL INVERS( XG(I, J) , YG(I, J) ,FX(I , J) ,FY(I, J) ,X, Y, 
1,EX,EY) 

WRITE (6,6) X, Y, EX, EY 

6 FORMAT (IH , 4 (17X, E13 . 5) ) 

GO TO 102 

103 CALL INVERS( XG (I , 1) , YG (I , 1) , 0 . , 0 . ,X, Y,EX, EY) 
WRITE(6,6) X,Y 
102 CONTINUE 

IF ( KK(I) .-LT. 0 ) PRINT 7 

7 FORMAT (/20X, 'SUBPROGRAN SCAN FAILED TO CONVERGE FOR 
1 -NEXT- POINT') 

101 CONTINUE 
STOP 
END 
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SUBROUTINE INVERS (U , V , EU , EV , X, Y , EX , EY) 
COMPLEX Z, W, S 
W = CMPLX(U,V) 

S - CSQRT(W**2 - 4) 

IF ( U ,LT. 0.0 ) S S 

P = REAL(S) 

Q = AIMAG(S) 

Z = (W + S)/2 

X = REAL(Z) 

Y = AIMAG(Z) 

G = ( X*P + Y*Q ) / ( P**2 + Q**2 ) 

H = ( Y*P - X*Q ) / ( P**2 + Q**2 ) 

EX = EU*G + EV+H 
EY = -EU*H + EV*G 
RETURN 
END 


The user must supply a data card to be read with the 
format (3I5,3F10.0) and containing the following; 

N; Number of data points not larger than 30, 

NPOTEN : Number of contours not larger than 10 , 

NPOINT: Number of points on each contour not 

larger than 50 , 

PMIN.; The first potential value among the 
contours to be calculated, 

PSTEP: The algebraic separation (potential- 

dimension) between contours, 

ADJ : Adjustment factor ranging approximately from 

0.1 to 1. The larger it is the greater the 
separation of two successive points on one 
contour becomes. 

The user must also supply N data cards (2F20.0) having the 
entries^ 

x’(I), V(I): Boundary coordinates and voltages 

on the straight line (u-axis in 
w-plane) . 
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3. Subroutine GREEN 


SUBROUTINE GREEN (M, NPOTEN, NPOINT, PMIN, PSTEP, ADJ) 
COMMON /ARR/ XG(10 , 50) , YG (10 , 50) , FX(10 , 50) ,FY(10, 

1 , 50),PP(10) ,KK(10) 

COMMON /ORI/ C(30), V(30) , BE(30) , AL(30) , PY, NN 
FY = 3,1415926 
DO 11 I = 1, NN . 

BE(I) = ( V(I+1) - V(I) ) / ( C(I+1) - C(I) ) 
AL(I) = V(I) - BE(I) * C(I) 

11 CONTINUE 
M = 0 

DO 12 I = 1, NPOTEN 
KEY = 0 

P = PMIN + (I-l) * PSTEP 
DO 13 J = 1, NN 
IF ( M .EQ. 10 ) RETURN 
ORDI =0.0 

IF ( (P-V(J))*(P-V(J+D) ) 131, 132, 13 

131 ABSS = (P-V(JX) ♦ ;(C(J+1)~C(J)) /'(V,(Jh-1)-V(J))+C(J) 
GO. TO 130 

132 IF( P .EQ. V(J) ) GO TO 133 
GO TO 13 

133 ABSS = C(J) 

130 M = M + 1 

KEY = 1 

SNS = - SIGN(ADJ,BE(J)) 

KK(M) = NPOINT 
PP(M) = P 

XG(M,1) = ABSS 
YG(M,1) = ORDI 

ORDI = ABS(C(J+1)-C(J)) / 200.0 
DO 14 K = 2, NPOINT 

CALL SCAN ( P, PSTEP, BE(J) , ABSS, ORDI, EX, EY, 

1 E2,KASE ) 

IF (KASE) 142,141,142 

141 XG(M,K) = ABSS 
YG(M,K) = ORDI 
FX(M,K) = EX 
FY(M,K) = EY 

ABSS = ABSS - SNS* ABS( PSTEP) *EY/E2 
ORDI = ORDI + SNS* ABS( PSTEP) *EX/E2 
IF(ORDI.LE.O.) KK(M) = K 

IF (ORDI) 13, 13, 14 

142 KK(M) = (K-1) * KASE 

GO TO 13 

14 CONTINUE 

13 CONTINUE 

IF ( KEY .EQ. 0 ) \YRLTE(6,7) P 
7 FORMAT (//20X, 'THERE IS NO CONTOUR CORRESPONDING TO 
1»P = ' ,E13.5//) 
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12 CONTINUE 
RETURN 
END 

SUBROUTINE SCAN (P,PSTEP,BE,X, Y, EX,EY,E2,KASE) 
DO 40 I = 1, 50 

EX = FIELDX (X,Y) 

EY = EIELDY (X,Y) . 

E2 = EX* *2 + EY**2 
TEST = BE**2 * l.OE-6 
IE ( E2 .LT. TEST ) GO TO 44 
DIE = POTEN (X,Y) - P 
CRI = ABS(PSTEP) * l.OE-4 
IF •( ABS(DIF) .LE. CRI ) GO TO 42 
X = X + DIF*EX/E2 
Y = Y + DIF*EY/E2 
IF(Y,LE.O.) GO TO 44 
40 CONTINUE 
KASE = -1 
RETURN 
42 KASE = 0 
RETURN 
44 KASE = 1 
RETURN 
END 


FUNCTION POTEN (X,Y) 

COMMON /ORI/ C(30), V(30) , BE(30) , AL(30) , PY, NN 
R(T) = (1/PY) * ATAN2(T-X,Y) 

S(T1,T2) = (0,5*Y/PY) * ALOG( (Y**2+ (T2-X) **2) 

1 /(Y**2+(T1-X)**2) ) 

N = NN + 1 

POTEN = 0.5 * ( V(l) + V(N) ) + V(l) * R(C(D) - 
1. V(N) * R(C(N)) 

DO 20 I = 1, NN 

20 POTEN = POTEN + ( BE(I)*X + AL(I) ) * ( R(C(I+1)) - 
1, R(C(D) ) + BE(I) * S(C(I) ,C(I+1)) 

RETURN 

END 


FUNCTION FIELDX (X,Y) 

COMMON /ORI/ C(30), V(30) , BE(30) , AL(30) , PY, NN 
R(T) = (1/PY) * ATAN2(T-X,Y) 

S1(T) = (Y/PY) / (T-X)**2 + Y**2 ) 

S2(T) = (T-X) * S1(T) 

N = NN + 1 

FIELDX = V(l) * S1(C(1)) - V(N) * S1(C(N)) 

DO 30 I = 1, NN 

30 FIELDX = FIELDX + ( BE(I)*X + AL(I) ) * ( S1(C(I+1)) 

1 - S1(C(D) ) + BE(I) * ( S2(C(I+1)) - S2(C(D) 

2 - R(C(I+1)) + R(C(D) ) 

RETURN 

END 
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5 / / 

FUNCTIONS FiELDY (X, Y) 

COMMON /OKI//C (30) , V(30) , BE(30) , AL(30) , PY, NN 
R1(T) = C (T-X)*+2 + Y**2 ) 

R2(T) = (T/X) / R1(T) 

S(T1,T2) ? 0.5 * ALOG( R1(T2) / Rl(Tl) ) 

N = NN^+/l 

FIELDY ? V(1) * R2CC(D) - V(N) * R2(C(N)) 

DO .35 / = 1, nN 

35 FIELDY = FIELDY + ( BE(I)*X + AL(I) ) * ( R2(C(I+1)) 

1 - R2(C(D) ) - BE(I) + ( Y**2 + (1/R1(C(I+D) 

2 -^/R1(C(I))) + S(C(I),C(I*1))'.) 

FIELDY = FIELDY / PY 

. RETURN 
END 


Subroutine GREEN is the implementation of Eqns-. (Ill- 

/ 

1) to (III-9) and the process described in Section B. This 


subroutine references three other subprograms that provide 
potential and field components at a particular point. 

After GREEN is executed, there are four two-dimensional 


arrays and two one-dimensional arrays in the block COMMON 
|ARR! available for printing or further operations. The 


argument M in the CALL statement indicates the -number of 


contours which GREEN has found where 0 < M < 10. GREEN 
"there is no contour corresponding to P" when it cannot be 
found. The array contents are. as follows: 


XG(.I,J) : 

til 

X coordinates of the J 
contour on the z plane. 

. . _th 
point, I 

YG(I,J) : 

■4“ 

y coordinates of the 

point , I^^ 


contour on the z plane, 

FX(I, J) : 

X component of the field 
YG(I,J)), 

at (XG(I,J), 

FY(I,J) ; 

y component of the field 
YG(I,J)), 

at (XG(I,J), 
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